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ABSTRACT 

Using a simplified model of cascade pair creation over pulsar polar caps presented in two previous 
papers, we investigate the expected gamma-ray output from pulsars' low altitude particle acceleration 
and pair creation regions. We divide pulsars into several categories, based on which mechanism truncates 
the particle acceleration off the polar cap, and give estimates for the expected luminosity of each category. 

We find that inverse Compton scattering above the pulsar polar cap provides the primary gamma 
rays which initiate the pair cascades in most pulsars. This reduces the expected 7-ray luminosity below 
previous estimates which assumed curvature gamma ray emission was the dominant initiator of pair 
creation in all pulsars. Even for the brightest pulsars where curvature radiation sets the height of 
the pair formation front (PFF), we find predicted luminosities too low to explain the EGRET pulsars, 
suggesting that the source of that emission is an outer magnetosphere accelerator. The predicted polar 
cap luminosities are large enough, however, to be observable by upcoming 7-ray instruments, which 
provides a firm test for this theory. 

Subject headings: Acceleration of particles — gamma rays: theory — pulsars: general 



1. INTRODUCTION 

Although only a few pulsars have so far been detected in 
7-rays, mostly via the EGRET telescope on CGRO (Ul- 
mer 1994; Thompson, et al. 1997), this number should 
steadily grow as a new generation of 7-ray instruments, 
such as INTEGRAL and GLAST, is brought online. This 
should greatly assist theoretical understanding of pulsars, 
since the highest energy photons provide a direct window 
into the underlying mechanisms thought to lead to pulsar 
emission of all types. 

The primary photon emission mechanisms important in 
radio pulsars arc inverse Compton scattering (ICS), in 
both the resonant (RICS) and Klein-Nishina nonresonant 
(NRICS) modes, and curvature emission. The relative im- 
portance of these mechanisms is still unclear. Most 7-ray 
pulsar papers have concentrated on the curvature emis- 
sion (Zhang & Harding 2000; Romani & Yadrioglu 1995; 
Daugherty & Harding 1982), while an increasing number 
of polar-cap physics papers have emphasized the impor- 
tance of ICS (Dermer 1990; Sturner 1995; Sturner, et al. 
1995; Luo 1996; Harding & Muslimov 1998). 

In the polar-cap acceleration model, particles are ex- 
tracted from the polar cap and accelerated by large 
rotation-induced electric fields, forming the primary beam. 
These particles then emit primary 7-ray photons due to 
ICS and curvature emission, and these photons interact 
with the pulsar magnetic field to create electron-positron 
pairs. The density of these secondary pairs increases with 
height as more and more photons pair-produce, until the 
pair density is sufficient to short out the accelerating field. 

Historically, this shorting of the accelerating field was 
thought to occur in a thin layer after the creation of the 
first pair; hence, the region of no pair creation, iSy ^ 0, 
was thought to be separated from the region of copious pair 
creation, i^y « 0, by a thin "pair formation front" (PFF). 



ICS photons, however, create small numbers of pairs at low 
altitudes, breaking the connection between the altitude of 
the first created pair and the altitude at which the elec- 
tric field disappears. Since the altitude of the first created 
pair has no inherent dynamical significance, we use the 
term PFF in this paper always as applying to the altitude 
at which the accelerating field is shorted. Other effects, 
such as the polarization of the generated pairs, will also 
begin at the point of first pair creation, but these effects 
must likewise reach some threshold before they affect the 
dynamics of the beam. 

Once the accelerating field has been shorted out, the 
primary beam coasts, continuing to emit 7-rays. The to- 
tal 7-ray output of a polar cap is then a combination of 
the synchrotron 7-rays produced by the created secondary 
particles, the ICS radiation emitted by the secondary par- 
ticles, and the primary 7-rays emitted by the primary 
beam. If the primary beam is radiation-reaction limited, 
this emission efficiently converts the beam energy into 7- 
rays; otherwise, only a small fraction is extracted. 

We find that the 7-ray output of pulsars falls into two 
categories. For the majority of pulsars, non-resonant ICS 
stops the beam acceleration at small Lorentz factors where 
7-ray emission is inefficient, leading to low luminosities. 
For the remaining pulsars, the beam is accelerated to high 
Lorentz factors, resulting in efficient 7-ray emission, and 
a high luminosity. 

Using the results of Hibschman & Arons (2001a) and Hi- 
bschman & Arons (2001b), henceforth Papers I and II, we 
examine the boundary between these categories of pulsars, 
and predict the luminosities and spectral characteristics of 
these objects. 

2. MODEL 

Inverse Compton scattering depends strongly on the 
temperature of the neutron star polar cap, making the 
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thermal cooling model eliosen for the neutron star an im- 
portant part of the theory. For simplicity, we will assume 
the temperature of the polar cap is entirely due to the 
polar cap heating discussed in Paper I. We use the accel- 
eration model of Muslimov & Tsygan (1992); Muslimov & 
Harding (1997), in the simplified form described in Paper 
I, and neglect spatial variations across the polar cap. 

This accelerating potential, expressed in units of mc^/e 
so as to give the expected particle Lorentz factor, is 



^high{t) = ^it 



mc^e-^, t<l 
mc'^e~^, t > 1 



(1) 



where $i = 5.14 x 10^Bi2P~^/^, t = s/si, and si = 
8.87 X lO^^P-i/^i?*, while .s is the altitude above the stel- 
lar surface. 

Although we use general relativity though the Muslimov 
and Tsygan accelerating potential, we neglect other rela- 
tivistic effects, such as the changes to the magnetic field 
near the surface and the gravitational red-shift of the emit- 
ted photons. These effects are of order 10-15%, so includ- 
ing them would pretend to greater accuracy than justi- 
fied. The accelerating potential is unique, as it is only the 
charge difference created by the rclativistic contribution 
which creates the starvation electric field. 

2.1. Emission rates 

First, we consider the emission of a single primary beam 
particle, moving along the field lines above the pulsar cap. 
We neglect, for now, discussion of the secondary particles 
created by pair-production, as the total energy emitted by 
these particles is clearly limited by that emitted by the pri- 
maries. For this section, we effectively assume that only a 
negligible fraction of the energy emitted in primary 7-rays 
remains in the generated electron-positron plasma. From 
the results of Paper II, this naturally follows in low-iJ pul- 
sars {B12 i 1) and results in other pulsars due to the RIGS 
of the secondary pair plasma. 

Since NRICS is only logarithmically-dependent on the 
Lorentz factor of the beam, the NRICS power emitted is 
limited primarily by the attenuation of the background 
thermal photons as the beam particles move away from 
the star. If we assume a hot polar cap of angular radius 
9c = y/ilR*/c, where f2 is the angular velocity of the pul- 
sar, O = 2'7r/P and ii* is the radius of the neutron star, 
assumed to be 10 km, the total energy radiated by one 
particle is 



Enr = 



OcR* 



Pnr = 6.1 X 10-4 p- l/2j^2 g^gg (2) 



The approximate form used for the NRICS power, Pnr, 
is that given in Paper I, and Tq is the temperature of the 
polar cap, in units of 10^ K. 

The power emitted via RICS, while only logarithmically 
sensitive to the decline in the thermal photon flux, de- 
creases quickly with increasing Lorentz factor. Because of 
this dependence, most of the RICS power is emitted at 
low altitudes where the Lorentz factor is still small. The 
expected total output per particle is 



-Pfl(7mm) = 4.5 X 10-5 P^/^Bi2T^^^ ergs (3) 



Efi = 



where "fmin = ^s/kT is the minimum Lorentz factor at 
which thermal photons are upscattered into resonance 



— 450 P'^/^T'g meters, assuming a 



with the beam, Smin is the altitude at which the beam par- 
ticles reach ^mini and Pr is the power emitted by RICS 
emission, as given in Paper I. In the acceleration model 

from Paper I, Smin 
star of radius 10 km. 

The power emitted by curvature radiation is strongly 
dependent on the Lorentz factor of the beam, varying as 
7^, so most of the energy emitted by curvature radiation is 
emitted as the beam coasts above the PFF. Once particle 
acceleration stops, the Lorentz factor of the primary beam 
declines according to 

7(S > SPFF) = IPFF 1 + 3 In — 

V IPFF C 1-1- SpFF J 

(4) 

assuming a dipolar magnetic field. The upper limit may 
be estimated by finding the altitude at which the high- 
energy primary beam decouples from the magnetic field. 
Equating the energy density in the beam to the energy 
density in the magnetic field gives a decoupling height of 

frnax = '^12'iB\l^ P^^^^j ^^'^P*, whcrc 77 is the Lorentz 
factor of the beam in units of 10^. The PFF is close to the 
surface, so the logarithm above is approximately 8.3. 
The total curvature energy emitted is then 

Ec = {iPFF - l{Smax))mC^ = 8.2(7pirF7 - n{Smax)) ergS. 

(5) 

Comparing the emitted energies, equations (2), (3), and 

(5), we find that for typical Lorentz factors of "Ipff > 10^, 
only the curvature emission may radiate any apprecia- 
ble fraction of the beam particle energy. The minimum 
Lorentz factor at which radiation reaction is important is 
then 

7ii^i = 2.33x 10^PV^ (6) 

which is the Lorentz factor at which the primary beam par- 
ticles lose half their energy to curvature radiation. Above 
this Lorentz factor, roughly all of the energy in the beam 
is lost to 7-rays; below, the beam propagates without sig- 
nificant radiation losses. 

The expected curvature energy loss in these two regimes 

is 



(7) 



-Pci^PFF) = 1.51 X 10-^^ p-^jj,pp ergs (8) 



Ec ~ JPFF rac? = 8.2 x IQ^'^^ppp ergs 
if radiation reaction is important and 

8.3P* 

Ec ~ - 

if not. 



2.2. Luminosity 

Using the polar cap model from Paper I, we can clas- 
sify pulsars according to the mechanism which sets the 
PFF and whether the beam is radiation reaction limited 
{iPFF > Irr)- The pair formation model then gives the 
altitude of the PFF, sppp, and the Lorentz factor at that 
altitude, jppp. 

Given jppp, we can compute the total expected lumi- 
nosity by multiplying the total energy emitted by a sin- 
gle beam particle by the number of particles emitted by 
the polar cap, N = ncjcOcRl = 1.37 x lO^oPiaP"^ s"\ 
where ncj is the expected Goldreich-Julian number den- 
sity, no J = ■ B/2-Kce. 



Gamma-ray emission 
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As a slight modification to the model in Paper I, wc find 
that, in comparison with the numerical results of Paper II, 
the PFF from curvature emission is more accurately found 
by finding the altitude at which the first pair is formed. 
This is due to the steadily increasing intensity of curvature 
emission with increasing Lorentz factor. 

A curvature photon emitted at altitude s will pair pro- 
duce at 



1 

s± = s + - 



4ec(7(s)) 



(9) 



where ec is the typical curvature photon energy, ec 



5.8 X 10 



-19^3 



and €a is the scaling energy for pair 



production from Paper I, Ca = 2166 B{^^P^/^/p. Here fp is 

the ratio of the actual field line radius of curvature to the 
radius of curvature of the dipole field line which intersects 
the stellar surface at 6c; for the remainder of the paper, 
this is taken to be 1. 

The minimum value of this is at 

spFF,c = 1.91 B^^P"^fl'^R., (10) 
if the PFF takes place in the linear regime of the acceler- 
ating potential at s > OcR*, or 



4/7pll/l4 



SPFF,C = 0.211 i?i7"P^^/^V;/'i?„ (11) 

if in the quadratic regime at s < 6'ci?*. 

Using the scmi-numcrical model of Paper II, we then 
classify pulsars by which the emission mechanism pro- 
duced the PFF, by whether the PFF occurred at low or 
high altitude (in comparison with the polar cap width), 
and by whether the beam was radiation-reaction limited. 
In principle, this gives 12 categories; in practice, there are 
only five important divisions. 

Out of the 540 pulsars in the Princeton pulsar cata- 
log with positive P, the majority, 315, had a PFF set 
by NRICS at high altitude with jpff < Irr- For all 
of these objects, curvature radiation is the primary en- 
ergy loss mechanism in the beam drift region above the 
PFF. Although most of the radiated energy is curvature 
emission, it is the relatively sparse, but individually much 
higher-energy, NRICS photons that first pair produce and 
set the PFF. 

In the model of Paper I, self-consistent cap heating pro- 
duces a PFF height and final Lorentz factor of 



S PFF. NR. 



pc 



lPFF,NR,pc ■ 



: 0.749 Bi"2'*/^P"/^i?* (12) 
: 4.34 X 10<5 B-//^P-V6 (13) 



which corresponds, through equation (8) to a total lumi- 
nosity of 

Lnr = 7.3 X 102S p-2'/-V-"/3 erg s'l. (14) 

This is much lower than previous estimates of pulsar 
7-ray luminosity, and is by far the most common case. 
However, this estimate is only strictly accurate for mag- 
netic fields less than or on the order of 10^^ Gauss due to 
the increasing fraction of particle energy left in the par- 
ticles at higher values of B, as discussed in Paper II. For 
higher fields, the luminosity will be larger, as NRICS is 
less efficient, but still less than the value expected from 
pure curvature radiation. 

The next most common case is high-altitude inertially 
limited RICS, accounting for 126 of the 540 pulsars. For 
most of these pulsars, the 7-ray emission in the beam drift 



region above the PFF is again dominated by curvature 
emission; the 23 of these where this is not true are among 
the least luminous pulsars, and so we neglect them. 

From Paper I, the PFF height and final Lorentz factor 
in this case are 

(15) 
(16) 



spFF,R,pc = 9.66 B^^*"'P^'^R, 



5.60 X 10^5-3^/^-3/2 



iPFFM.pc 

which corresponds, via equation (8), to a total luminosity 
of 

L^^^^ = 2.0 X 10^3 B-^^I^P-^ ergs s'^. (17) 

The third category combines those pulsars where either 
curvature or NRICS sets the PFF at high altitude and 
where the beam is radiation reaction limited, for a total of 
57 objects. In practice, we find that these pulsars are all 
well-modeled by the expected PFF for curvature emission, 
equation (10), and equation (7) for the energy loss. This 
gives a final Lorentz factor and total luminosity of 



7^^;;c = l-10x lO^p-i/4 
L^'s'* = 1.2 X 1031 P12P-3/4 ergs s"! 



(18) 
(19) 

which is essentially the same as the Zhang & Harding 

(2000) value for their regime II. As they mentioned, this 

1 /2 

preserves the empirical L-y oc Lg^, relation. These are 
among the brightest 7-ray pulsars, behind only category 4 
in total luminosity. 

The fourth category is similar to the third, in that it 
includes radiation reaction limited beams where the PFF 
is set by both curvature and NRICS emission, but at low 
altitudes rather than high. These are the 21 brightest, 
highest-potential pulsars. Using equations (11) and (7) 
gives a Lorentz factor and luminosity of 



IPFFfi = 2-90 X 10^Pi2 
L 



14 



3.3 X 1031 P;^/ V-27/14 g^gg g-1 



^PFF,R,pc ■ 

,,low 



:4.03Pf2^p3/^P* 



rPFF,R,pc^^-OQxw'°B^ip-^/\ 



(20) 

— ^-^ ^ -"12 ^ " (21) 

which corresponds to regime I of Zhang & Harding (2000) . 

The final category consists of the 13 pulsars where RICS 
sets the PFF in the low-altitude regime, with an inertially 
limited beam. These are all pulsars with fields well in ex- 
cess of 10^3 Gauss, field strengths beyond the expected 
range of validity of this theory. The the low-altitude PFF 
height and Lorentz factor are 

(22) 
(23) 

For these pulsars, RICS emission produced more gamma- 
rays than curvature, so using the emitted energy from 
equation (3), we find a luminosity of 

L'^^ = 7.2 X 10^6 P^/^p-67/32 g^gg g-1 (24) 

As in Paper I, the emission mechanism with the small- 
est PFF height is the dominant mechanism. Owing to the 

multiplicity of categories, generalizations are difficult, but 
we can derive a few formulae for the boundaries between 
regimes by considering only the high-altitude limits. 

In this simplest appoximation, the boundary between 
the NRICS-dominated and curvature-dominated pulsars 
is at P12 = 0.061 pi/**. Pulsars with weaker magnetic 
fields are dominated by curvature, and equation (19) is the 
appropriate luminosity, while pulsars with fields stronger 
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than this are dominated by NRICS and equation (14) is ap- 
propriate. However, this ignores the weakening of NRICS 
with increasing magnetic field discussed in Paper II, and 
so should be taken as a statement that the millisecond 
pulsars arc certainly controlled by curvature, while the 
higher-field, longer-period pulsars are favored by NRICS, 
but must be examined carefully, using either the full alge- 
braic results from Paper I or the semi-numerical model of 
Paper II. 

The boundary between NRICS and RICS lies at 

— 5/7 

Pnr,r = 6.81 i?i2 5 with higher periods favoring RICS, 
while that between RICS and curvature is at Pr,c = 
again with RICS dominating at longer pe- 
riod. For pulsars with magnetic fields larger than approxi- 
mately 4 X 10^^ Gauss, NRICS is ineffective in setting the 
PFF, and the only active mechanisms arc RICS and Cur- 
vature; at mid-range fields between 3 x 10^^ Gauss and 
4 X 10^2 Gauss, NRICS sets the PFF, while at lower fields, 
curvature sets it. 

2.3. Flux 

Since in this model the observed 7-rays originate from 
particles moving along the field lines above the pulsar po- 
lar cap, the predicted luminosities translate directly into 
a predicted flux. Since the polar cap is the source of 
the 7-rays, they are beamed into a cone of opening angle 
9 » (3/2)^c- In general, this opening angle will vary with 
the altitude of the emission, but for most pulsars the al- 
titude of the PFF is small compared to the stellar radius, 
so simply evaluating the opening angle at the surface is 
sufficient. If the luminosity is spread evenly throughout 
this cone, this produces a peak flux of 

^p^^^ = ^ (25) 
where d is the distance to the pulsar. This is the flux seen 
while looking "down the barrel of the gun." Only if the 
spin axis of the pulsar is nearly aligned with the line of 
sight will the observed flux be on this order. The aver- 
age flux is reduced by O/tt if the spin axis is perpendicular 
to the line of sight and by approximately ^/tt sin a in the 
general case, where a is the angle between the magnetic 
moment and the rotation axis and we have effectively as- 
sumed that the observer's line of sight passes through the 
center of the emitting cone (/3 = 0). This represents the 
fraction of the pulse period where the emitting cone is di- 
rected towards the observer. 

= ^ 20 ■ 32- (26) 

This average flux is the observable, not the total lu- 
minosity. In terms of the luminosities from the previous 
section, this flux is 

<i)ave = 3.7 X 10-^^P^/^{sma)-^d^^^Levgs cm-2 s-\ 

(27) 

Since the sin a term is a factor of order unity, and these 
approximations are only good to approximately a factor 
of 2, we simply set it to 1 for the remainder of the pa- 
per. We plot the expected observable flux as a function 
of the pulsar period using a fiducial pulsar distance of 1 
kpc in Figure 1, to show the inherent dependencies on 
pulsar parameters, and using the estimated distances of 
each individual object in Figure 2, to show the predicted 
observable flux. 



3. SPECTRAL SHAPE 

In the previous section, we discussed the total energy 
output of the pulsar; here, we turn to the expected shape 
of the spectrum itself. In Paper II, we found that a cascade 
of pair creation from a single absorbed photon produces 
a response with a power law index of —3/2, due to the 
reprocessing of synchrotron photons. In order to be ab- 
sorbed, however, the photon must have an energy greater 
than €min = 5134i?{^^P^/^mc^, presuming that the pho- 
ton was emitted parallel to the magnetic field at the edge 
of the polar cap at the surface of the star. 

Since most of the power emitted in 7-rays from pulsar 
polar caps is due to curvature emission, except for the few 
extreme high-field objects where RICS dominates, the en- 
ergetics of the curvature photons determine the shape of 
the spectrum. The minimum Lorentz factor for curvature 
pair production may be found by equating ec(7) = emini 
yielding 

7c = 2.01 X 10^ B~2^^P^/^ (28) 



which is the Lorentz factor at which the critical curvature 
energy ec(7) = 5.8 x 10~^^P8 ^7'^rnc^ equals the minimum 
energy to pair produce, Cmm, assuming a dipole field ra- 
dius of curvature. 

If ^pFF > 7c : then the copious curvature photons will 
pair-produce, and the observed radiation will be the syn- 
chrotron emission of the generated pairs. At low energies, 
e < €0/(1 + a^) In A in the notation of Paper II, the spec- 
tral index is the —2/3 of unprocessed synchrotron radi- 
ation, while at higher energies, the spectral index is the 
characteristic —3/2. 

If 7pff < 7c, then the curvature photons will be ob- 
served directly. The beam energy loss in this case, for all 
pulsars examined, is small enough that the characteristic 
curvature energy remains effectively unchanged. There- 
fore, the unmodified —2/3 spectrum of curvature radia- 
tion will be seen, extending from low energies up to ecij)- 
This energy will clearly be lower than the traditional es- 
timate of the cut-off energy of and dependant on the 
mechanism which sets the PFF. 

Figure 3 shows the results of dividing the pulsars into 
two categories, based on this division. Out of the top 
20 brightest expected polar-cap 7-ray pulsars, only two 
are expected to have the —2/3 power law, J0953-I-0755 
and .11932-1-1059. To illustrate the difference, we plot 
in Figure 4 the numerically-calculated 7-ray spectrum of 
J1952-I-3252, which is expected to have a —3/2 spectrum, 
and J1932-I-1059. The numerical method used was that 
described in Paper II, and the results confirm the expec- 
tations. 

The maximum energies of these spectra depend on 
whether the beam is radiation-reaction limited or not, not 
on the minimum energy for pair production, Ca . All photon 
energies higher than Ca are converted by pair production 
to photons of lower energy, but steadily increases with 
altitude as r^, eventually surpassing the maximum energy 
of the raw photon spectrum emitted by the beam parti- 
cles. At that point, and beyond, the maximum energy of 
the spectrum is set by the energy of the beam itself. 
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Fig. 1. — 7-ray flux expected from pulsars, as a function of pulse period, if all pulsars are placed at a fiducial distance of 1 kpc. 
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Fig. 2. — 7-ray flux expected from pulsars, as a function of pulse period, given the actual pulsar distances. 
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Spectral classes 
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Fig. 3. — Expected observed 7-ray flux vs. pulsar period. Each pulsar has been classified into those where curvature radiation is expected 
to produce pairs {crosses) and those where it is not (circles). Pulsars where curvature emission produces pairs are expected to have a spectral 
index of —3/2, while those which do not are expected to have an index of —2/3. 



10^ 



PSR1 952+3252 




Fig. 4. — Predicted 7-ray spectrum produced by a single primary particle for 1952+3252, illustrating the strong curvature-induced pair 
production regime, with the resultant u = —3/2 power-law spectrum. The spectrum is given in units of N jmr? . 
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Fig. 5. — Predicted 7-ray spectrum produced by a single primary particle for 1932+1059, illustrating the sparse pair production regime, 
with the resultant v = —2/3 power-law. The spectrum is given in units of N/mc^. 



Since curvature radiation is the strongest emission 

mechanism for ah of the brightest pulsars and all but a 
small minority of the others, the maximum energy will be 
the characteristic curvature energy, evaluated at the high- 
altitude coasting Lorentz factor of the beam. However, 
the relevant Lorentz factor is slightly different from that 
calculated in the previous section, equation (6). Since the 
radius of curvature of the field steadily increases with alti- 
tude, the curvature energy of a coasting beam will steadily 
decrease, as r~^/^. The maximum energy observed will 
then arise from emission from an intermediate regime; 
high enough that the magnetic field no longer absorbs pho- 
tons, but low enough that the radius of curvature is still 
small. In this case, the appropriate limiting Lorentz fac- 
tor is where the scale height for energy loss is equal to the 
stellar radius, 7^^ = 3.56 x Iff P^/^. 

If the beam is radiation-reaction limited in this sense, 
i.e. if 7pff > i'rri then the maximum energy is the cur- 
vature energy evaluated at ^'nn, or 

e^ax = ec(7kfl) = 10.2P^/'GeV (29) 

where we have evaluated the curvature energy at a radius 
of 2i?*. 

If the beam is not radiation-reaction limited, then the 
maximum energy depends on the Lorentz factor of the 

beam, through 

= ecilPFF) = 3.2 X 10-22 p-"^^lpp GeV. (30) 

The expected observable fluxes and maximum energies 
for the 25 brightest pulsars are shown in Table 1. 

4. NUMERICAL MODEL 

Using the full numerical system described in Section 4 of 
Paper II, we have run several simulations of these objects. 
These results confirm our conclusions about the diff'erent 



regimes of pair production discussed above, with the nu- 
merically calculated luminosity remaining within approx- 
imately a factor of 3 of the simple analytic model, with 
the numerical model always lower, due to the effects of 
pair creation on the tails of the distribution. Over the 25 
brightest pulsars, the numerical results are on average of 
a factor of 2.1 lower than the theoretical results, although 
the brightest pulsars show more variance. 

The calculated pulsar spectra matched expectations, al- 
though in several cases a low-intensity high-energy tail 
due to NRICS was observed, extending up to -fPFF with a 
power law exponent of roughly -2. In general, this tail con- 
tains a negligible portion of the energy, namely the energy 
predicted by equation (2). 

We also ran the numerical model using alternative heat- 
ing models. We found that, since in the brightest objects 
curvature emission sets the PFF, the precise temperature 
model mattered little for the observed objects. Using dif- 
ferent models of the stellar temperature only changed the 
expected luminosity by on the order of 20%. However, for 
the lower-luminosity ICS-dominated pulsars, the temper- 
ature is far more important, as the ICS process relies on 
the thermal photon bath, potentially allowing future 7-ray 
observations to discriminate between thermal models. 

5. DISCUSSION 

In Table 2, we compare the observed flux, the flux by 
the model of Zhang & Harding (2000), the flux predicted 
by the semi- numerical model of this paper, and the flux 
computed by running the full numerical cascade model for 
each of the observed 7-ray pulsars. The flux predicted for 
both models was derived from the total luminosity using 
equation (27). 

The model of Zhang & Harding (2000) assumed that 
a = 30, decreasing the expected flux, while invoking a 
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Table 1 
Brightest pulsars 





P 




B 


r r 


ixrv 


Curv. 




Predicted Flux 






(Gauss) 


inGctl. 


limitGd? 


pairs? 


(_ijev j 


/ _2 

(^erg cm 


s ) 


J0835-4510 


0.089 


6.8 


X 


1012 


Curv 


Yes 


Yes 


3.07 


1.6 X 10- 


-8 


J0633+1746 


0.237 


3.3 


X 


1012 


NRICS 


Yes 


Yes 


3.73 


9.5 X 10" 


-9 


J0437-4715 


0.006 


1.2 


X 


109 


Curv 


Yes 


Yes 


0.78 


5.7 X 10- 


-9 


J0534+2200 


0.033 


7.6 


X 


1012 


Curv 


Yes 


Yes 


1.88 


4.8 X 10" 


-9 


J1932+1059 


0.227 


1.0 


X 


1012 


NRICS 


No 


No 


0.81 


1.1 X 10" 


-9 


J 1709 -4428 


0.102 


6.3 


X 


1012 


Curv 


Yes 


Yes 


3.29 


9.2 X 10" 


-10 


J0953+0755 


0.253 


4.9 


X 


lOii 


NRICS 


No 


No 


0.55 


6.2 X 10^ 


-10 


Ji3UU + i/4U 


U.OOd 


1.7 


X 


109 


Curv 


Yes 


Yes 


O.oi 


3.7 X 10- 


-10 


J1952+3252 


0.040 


9.7 


X 


lOii 


NRICS 


Yes 


Yes 


2.04 


3.3 X 10- 


-10 


J1048-5832 


0.124 


7.0 


X 


1012 


Curv 


Yes 


Yes 


3.61 


2.9 X 10 


-10 


J1012+5307 


0.005 


5.6 


X 


108 


Curv 


Yes 


Yes 


0.74 


2.3 X 10" 


-10 


J2043+2740 


0.096 


7.0 


X 


lOii 


NRICS 


Yes 


Yes 


3.18 


2.0 X 10 


-10 


J0742-2822 


0.167 


3.4 


X 


1012 


Curv 


Yes 


Yes 


4.19 


1.9 X 10 


-10 


J0034-0534 


0.002 


2.3 


X 


108 


Curv 


Yes 


Yes 


0.44 


1.9 X 10" 


-10 


J1826-1334 


0.101 


5.6 


X 


1012 


Curv 


Yes 


Yes 


3.27 


1.7 X 10" 


-10 


J1803-2137 


0.134 


8.6 


X 


1012 


Curv 


Yes 


Yes 


3.75 


1.6 X 10" 


-10 


J1730-2304 


0.008 


7.9 


X 


108 


Curv 


Yes 


Yes 


0.93 


1.6 X 10" 


-10 


J0117+5914 


0.101 


1.6 


X 


1012 


NRICS 


Yes 


Yes 


3.27 


1.4 X 10" 


-10 


J1959+2048 


0.002 


3.3 


X 


108 


Curv 


Yes 


Yes 


0.41 


1.4 X 10- 


-10 


J1801-2451 


0.125 


8.1 


X 


1012 


Curv 


Yes 


Yes 


3.63 


1.3 X 10' 


-10 



Table 2 

Observed and predicted pulsar fluxes 





B 


PFF 






4'zhang 


b 


(ppred 




(pnurn 




Pulsar 


(Gauss) 


mech. 


(erg cm^2 


s-i) 


(erg cni^2 


s-i) 


(erg cm^2 


s-i) 


(erg cm^2 


s-i) 


J0534-h2200 


7.6 X 1012 


Curv 


1.3 X 10" 


-8 


5.9 X 10- 


-9 


4.8 X 10" 


-9 


1.8 X 10- 


-9 


J0835-4510 


6.8 X 1012 


Curv 


9.9 X 10- 


-9 


2.3 X 10" 


-8 


1.6 X 10" 


-8 


5.0 X 10- 


-9 


J0633+1746 


3.3 X 1012 


NRICS 


3.9 X 10- 


-9 


2.6 X 10" 


-8 


9.5 X 10- 


-9 


4.5 X 10- 


-9 


J1709-4428 


6.3 X 1012 


Curv 


1.3 X 10- 


-9 


1.3 X 10" 


-9 


9.2 X 10- 


-10 


3.0 X 10- 


-10 


J1513-5908 


3.1 X 1013 


RICS 


8.8 X 10- 


-10 


5.3 X 10- 


-10 


4.0 X 10' 


-13 


2.0 X 10- 


-13 


J1952+3252 


9.7 X IQii 


NRICS 


4.3 X 10- 


-10 


5.1 X 10- 


-10 


3.3 X 10" 


-10 


1.8 X 10- 


-10 


J1057-5226 


2.2 X 1012 


NRICS 


2.9 X 10- 


-10 


2.5 X 10- 


-10 


7.6 X 10- 


-11 


5.1 X 10- 


-11 


J1048-5832 


7.0 X 1012 


Curv 


2.5 X 10- 


-10 


4.2 X 10- 


-10 


2.9 X 10- 


-10 


7.5 X 10- 


-11 



'"Thompson, et al. (1999) 
''Zhang & Harding (2000) 



Gamma-ray emission 
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proposed ICS instability to move the acceleration region 
off the surface of the star. These two effects roughly can- 
cel, leaving the predictions of their model comparable to 
those of this model, which effectively uses cos a = 1 and 
acceleration near the surface. Due to the small numbers of 
reversed particles, we find no instability in the acceleration 
zone, especially in the cases where, in the language of this 
paper, curvature radiation sets the PFF. Hence we find no 
reason to raise the altitude of the acceleration zone. 

The numerically calculated result is substantially 
smaller than either the semi-numerical result or that of 
Zhang & Harding (2000), due to two major effects. First, 
the numerically-calculated PFF is lower than that pre- 
dicted by the analytic model by approximately 25% on 
average, which, since these pulsars operate in the low-level 
quadratic portion of the accelerating potential, reduces the 
Lorentz factor of the beam to approximately 60% of its an- 
alytically calculated value. This occurs because photons 
on the exponential tail of the curvature spectrum pair pro- 
duce and create a sufficiently dense plasma to short out the 
accelerating electric field at lower altitudes than expected 
in cruder calculations (Arons & Scharlemann 1979). 

The second reason is that some of the primary photon 
energy remains in the generated pairs, rather than being 
re-radiated. Both the analytic model of this paper and 
that of Zhang & Harding (2000) assume that all of the en- 
ergy emitted by the primary beam is eventually re-emitted 
in 7-rays, due to the combination of synchrotron emission 
from the created pairs and RICS extracting any remaining 
energy. In the numerical model, however, we calculate the 
pair spectrum itself and can determine what fraction of 
the energy in the pairs is re-radiated by RICS. 

As a quick approximation to the effects of RICS on the 
pair spectrum, we assume that all particles with an energy 
loss length scale equal to or less than the stellar radius 
re-emit all of their energy as lower-energy 7-rays. The 
minimum Lorentz factor is the lowest at which thermal 
photons could be scattered into resonance with the field, 
7mm = es/A/xfcT = 134.5 Bi2Tg"^A/x"^ The maximum 
Lorentz factor is the point where the particle's energy loss 

Q 1/2 

scale is a stellar radius, jmax = 4.05 x 10 -Bi2T'g . 

With these corrections, we see that the predictions of 
the polar cap model are low compared to the observations 
for most of the observed 7-ray pulsars. The difference is 
large for the Crab, for which the prediction is low by a fac- 
tor of 10, and for the high-field object 1513-5908, which is 
low by three orders of magnitude. The other objects are 
typically low by a factor of three. 

The high-field object 1513-5908 clearly deserves further 
examination. Not only is this magnetic field of this pulsar 
so large that the applicability of this model is question- 
able, due to high-field effects as discussed by Harding, et 
al. (1997), but preliminary studies of the spatial variation 
of 7-ray emission across the polar cap suggest that the 
core of this pulsar's beam should be much brighter than 
the edge field line considered in this model. 

The neglect of spatial variation across the polar cap lim- 
its the accuracy of these results; simple estimates suggest 
that including those effects would increase the expected 
flux by 50%, but more detailed study is required to be 
more concrete. Physically, the general-relativistic acceler- 
ation is strongest at the center of the pulsar, while the field 



line radius of curvature is smallest at the edge. Due to the 
smaller radius of curvature and gentler acceleration, pair 
production due to ICS processes is far easier at the edges 
of the polar cap than at the center, while pair production 
via curvature radiation is more likely on the central field 
lines. Together, these effects combine to place the pair for- 
mation front at a higher altitude in the center of the polar 
cap and lower near the edges, so that the central field lines 
are comparatively brighter than the edges. 

These effects should help raise the predicted luminosi- 
ties closer to those observed, bringing them within the 
expected margins of error for this study. Further work on 
the variation of the 7-ray emission across the polar cap 
is certainly required before any firm conclusions can be 
drawn. 

However, the polar cap model gives specific predictions 
for the beaming which are not consistent with the obser- 
vations. First, the 7-ray peaks should be in phase with the 
radio emission, which is not typically the case. Second, for 
all of the observed pulsars, the curvature emission energy 
loss, equation (4), predicts that 40-60% of the emitted en- 
ergy for these pulsars should be emitted within the first 10 
km above the surface. This energy would go into a cone of 
angular width eem^t = 2(3/2)2^/^0, or 0,^u = 3.5 P-^/^ 
degrees. This is far more focused than the observed pul- 
sars, which have double beam profiles spread over rota- 
tion phase 0.4, typically. If the observed pulsars were all 
oriented with the magnetic axis, rotation axis, and line of 
sight all roughly parallel, such broad profiles could be gen- 
erated by a polar cap model. Specific modeling of the polar 
cap model's beaming properties require the angle between 
the rotation axis and the magnetic axis to be less than 
45 degrees and the intrinsic opening angle of the beam to 
be as large as 30 degrees (Harding & Daugherty 1998). 
The restriction on obliquity is a priori improbable, and 
is known to cause difficulty with the population statistics 
of pulsars (Romani 1996), while the large opening angle 
for the intrinsic emission beam requires adding new fea- 
tures to the basic dynamics of the polar cap model not 
supported by our calculations. 

Given these beaming problems with the polar cap 
model, the outer gap models of Romani & Yadrioglu 
(1995); Yadrioglu & Romani (1995); Zhang & Cheng 
(1997) remain strong contenders for the EGRET pulsars. 
These models easily produce beam profiles which resemble 
the data, but are energetically more difficult. For exam- 
ple, recent work by Hirotani & Shibata (2001) addresses 
the over-abundance in TeV emission, but the emissivity in 
the GeV range remains an open question. 

However, even if the observed 7-ray emission originates 
in the outer magnetosphere, the flux predicted to arise 
from the polar cap itself is not so far less than the observed 
fluxes as to be unobservable in the near future. With the 
advent of the new 7-ray instruments, the threshold of sen- 
sitivity should be low enough that the polar cap 7-rays can 
be observed. Detection of these gamma rays, seen when 
looking down the barrel of the radio gun, and comparison 
of their emission phase and profile with the radio coun- 
terpart, would be a firm test of the long established, but 
not directly tested, theory of polar cap pair creation which 
underlies 30 years of theorizing about the origin of pulsar 
radio emission. 
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6. CONCLUSION 

The primary conclusion of this paper is that the cur- 
rently observed 7-rays pulsars are not representative of the 
bulk of pulsars. Currently, only the brightest 7-ray pulsars 
are observed; these pulsars are selected to be those where 
curvature radiation operates efficiently, which is not true 
of pulsars in general. 

Most pulsars have a beam energy controlled by nonrcs- 
onant inverse Compton scattering (NRICS), rather than 
curvature emission, and will emit significantly fewer 7- 
rays. Locating examples of these objects will be a chal- 
lenge to the new generation of 7-ray telescopes. The 
fainter the objects which can be seen, the stronger the dif- 
ferences should be between the NRICS cascades and the 
predictions of the curvature-based PFF models. 

In these fainter objects, the separation of the mecha- 
nism producing the PFF (and thus shorting out the accel- 
erating potential) from the mechanism producing most of 
the observable 7-rays should reveal itself through differing 
spectral indices. Due to the effects of inverse Compton 
scattering, there should exist two classes of observed 7- 



ray emission, one a raw curvature spectrum characterized 
by a spectral index of —2/3, and the other a saturated 
synchrotron response with a spectral index of —3/2. 

For the currently observed 7-ray pulsars, the predicted 
fluxes are consistently low, although, due to the neglect of 
spatial variation across the polar cap, they are within the 
accuracy of this study. The beaming issues with the polar- 
cap model remain; for this reason, the outer-gap models 
are still strong contenders for the observed pulsars, despite 
the energetic problems of those models. Even if the cur- 
rently observed emission is from an outer gap, improved 
7-ray sensitivity should reveal the signal of the polar caps, 
which is at worst only a factor of 5 below the observed 
fluxes. 

Further observations should also reveal the differences 
in the luminosity classes of polar cap 7-rays, according 
to which mechanism which creates the PFF, and the dif- 
ferences between the two spectral regimes, providing a 
straightforward test of this model. 

The author would like to thank the referee for several 
helpful suggestions. 
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